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The photonuclear contribution to charged lepton energy loss has been 
re-evaluated taking into account HERA results on real and virtual photon 
interactions with nucleons. With large Q 2 processes incorporated, the average 
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muon range in rock for muon energies of 10 9 GeV is reduced by only 5% 
compared with the standard treatment. We have calculated the tau energy 
loss for energies up to 10 9 GeV taking into consideration the decay of tau. 
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I/"") ■ A Monte Carlo evaluation of tau survival probability and range show that at 

energies below 10 7 — 10 8 GeV, depending on the material, only tau decays are 
important. At higher energies the tau energy losses are significant, reducing 
the survival probability of the tau. We show that the average range for tau is 
' shorter than its decay length and reduces to 17 km in water for an incident 

tau energy of 10 GeV, as compared with its decay length of 49 km at that 
^ ■ energy. In iron, the average tau range is 4.7 km for the same incident energy. 



I. INTRODUCTION 

Neutrino telescopes have the potential for detecting distant sources of high energy neutri- 
nos, for example, from Active Galactic Nuclei and Gamma Ray Bursters M. Upward-going 
muons from muon neutrino conversions to muons via charged current interactions with nu- 
clei are the main signal of muon neutrinos @,[|. In addition, muons as well as neutrinos 
are produced in the atmosphere by cosmic ray interactions with air nuclei. Underground 
detector measurements of muon fluxes as a function of zenith angle are one way to determine 
the atmospheric muon flux as a function of energy. At high energies, one expects the muon 
flux to reflect the onset of contributions from the production of charm in the atmosphere [|J . 
A good understanding of the muon energy loss at high energies is an essential ingredient to 
the analysis of the high energy atmospheric muon flux. 

Recent SuperK measurements of atmospheric neutrinos strongly suggest — > v T oscil- 
lation H with large mixing angle of sin 2 29 > 0.84 and a neutrino mass squared difference 
of 2 x 1CT 3 eV 2 < Am 2 < 6 x 10~ 3 eV 2 ||. Assuming z/ M — > v T oscillation of extragalactic 
neutrinos with SuperK parameters, about half of the muon neutrinos are converted to tau 
neutrinos on the way to the Earth. This leads to comparable fluxes of ultrahigh energy tau 
neutrinos and muon neutrinos at the Earth. Tau energy loss affects tau neutrino propagation 
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in the Earth, where tau neutrinos interact with nucleons to produce taus which subsequently 
decay [/]]. An understanding of tau energy loss at very high energies could help with the 
interpretation of long tracks produced by charged particles in large underground detectors. 

Future neutrino telescopes such as AMANDA, NESTOR, ANTARES and ICECUBE are 
aimed at detecting high energy events from extragalactic neutrino sources. The high energy 
behavior of muon and tau interactions with water or rock nuclei have implications for event 
rates and the eventual unfolding of their respective parent neutrino flux. 

One way to characterize the energy loss of charged leptons is to consider the average 
lepton energy loss per distance traveled (X in units of g/cm 2 ), which can be expressed in 
the form 

-(§) = a + P E (1.1) 

where E is the lepton energy, a is nearly constant, determined by the ionization energy 
loss, and (3 = J2i P l is weakly dependent on energy and due to radiative energy loss through 
bremsstrahlung, pair production and photonuclear scattering. Generically, 

N rymax daHv E) 

P\E) = tL d y y^l^l (L2 ) 

where y is the fraction of lepton energy loss in the radiative interaction: 

(1-3) 

for final lepton energy E'. The superscript % denotes bremsstrahlung (brem), pair production 
(pair) and photonuclear (nuc) processes. Avogadro's number is N and the atomic mass 
number of the target nucleus is A. For rock, typical values for initial muons of energy ~ 100 
GeV are a ~ 2.4 MeV/(g cm" 2 ) and (3 ~ 3 x 10- 6 /(g cm" 2 ) §. 

At low energies, continuous energy loss by ionization dominates muon propagation, but 
at higher energies (above ~ 10 3 GeV), losses through pair production, bremsstrahlung and 
photonuclear interactions dominate. In the case of the muon, pair production is the most 
important mechanism, but for taus, the photonuclear process is at least as important as pair 
production. The high energy extrapolation of the photonuclear cross section has the largest 
theoretical uncertainty in the contributions to energy loss. 

The experimental interest in high energy muons and taus make it useful to reanalyze 
the high energy photonuclear differential cross section. Bottai and Perrone 0] have reeval- 
uated the photonuclear contributions using the HERA results for a(jN). In this paper, we 
evaluate the photonuclear differential cross section using HERA results for real and virtual 
photon-nucleon scattering. Here, we use a parameterization of the measured electromagnetic 
structure function F 2 which agrees well with data over the full range of Q 2 : from real photon 
interactions with nucleons to highly virtual photon interactions with nucleons. One expects 
that the inclusion of the Q 2 ^> part of the photonuclear cross section will contribute to 
larger muon energy losses at high energy. Indeed, this is what we find at the highest muon 
energies considered. We use the same expression for photonuclear interactions of taus with 
the correction for the different lepton mass. 

In the next section and Appendix A, we outline the formulas for evaluating charged 
lepton energy loss, and we review the standard treatment of the photonuclear cross section 
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discussed by Bezrukov and Bugaev in Ref. [[K]]. We outline the Monte Carlo program for 



energy loss (and decay, for the tau). In the following section, with details in Appendix B, 
we describe our treatment of the photonuclear cross section using F 2 . In Section IV, we 
evaluate and compare our energy losses with the standard treatment in the case of muons. 
Section V includes f3 % , range and survival probability for taus taking into consideration the 
decay probability. In Section VI, we summarize our results and discuss some implications 
of the new high energy photonuclear contributions to charged lepton energy loss. 



II. CHARGED LEPTON ENERGY LOSS 

A standard compilation of the energy loss parameters for muon energy loss was per- 
formed by Lohmann, Kopp and Voss (LKV) ||. The standard formulas are reproduced in 
Appendix A. In evaluating the energy loss due to ionization, the Bethe-Bloch formula [TT| 



is used with parameters listed in Table I [12]. Bremsstrahlung energy losses are computed 



via the differential cross section of Petrukhin and Shestakov ||13||. Improvements of the 



bremsstrahlung calculation |14[ affect the muon intensity-depth relationship by only a few 
percent [[Uj and should affect the tau energy loss even less. The pair production differential 
cross section is parameterized by Kokoulin and Petrukhin |IE| . 

The standard photonuclear differential distribution used in muon energy loss calculations 
is by Bezrukov and Bugaev (BB) flIU| . Their differential cross section is based on a gener- 
alized vector dominance model with off-diagonal contributions, and the nuclear shadowing 
is evaluated using an optical model. The BB parameterization of the photonuclear cross 
section is a function of A^-^a^^^E), as exhibited in Eqs. (A10) and (All). Our alternative 
to the BB differential cross section is described in detail in the next section. 

For quantitative comparisons of the two approaches to photonuclear energy loss, we use 
a one-dimensional Monte Carlo propagation program to evaluate survival probabilities of 
muons and taus as a function of energy and depth X. By using the dE/dX formula and 
integrating, one can find the "range of average energy loss" R/ae)- Lipari and Stanev (LS) 
have shown that this is different than the average range (R(E)) |]17|| . With a stochastic 
treatment of bremsstrahlung, pair production and photonuclear interactions, fluctuations 
have the effect of decreasing the average range relative to the range of average energy loss. 
Accounting for fluctuations has implications for the downward muon rates where the tails 
of the survival probabilities are important |17[| . The Lipari-Stanev evaluation of the muon 



range relies on the Bezrukov-Bugaev photonuclear contribution. 

The one- dimensional approximation, where the outgoing charged lepton travels in the 
same direction as the incident lepton, should be adequate for the high energies considered 
here. More elaborate (three-dimensional) evaluations have been performed by other authors, 
for example, in the program MUSIC by Antonioli et al. ||15|| . The MUSIC program compares 
well with the one-dimensional Lipari-Stanev survival probabilities, both of which use the 
differential cross sections outlined in LKV. 

In constructing our Monte Carlo program, we follow the standard procedure of splitting 
the radiative energy loss into two terms, a continuous "soft" term for y < y cut , and a 
stochastic "hard" term for y cut < y < 1 in 

dE - a+ N E ( VcUt dyy da + N E ( l dyy d ° (2 1) 

dX A Jo dy A J ycut dy ' 
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while the ionization is treated continuously. We have chosen y cut = 10 3 following Ref. 
For taus, we have added a stochastic term for tau decay. 
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III. PHOTONUCLEAR CROSS SECTION USING F 2 

The photonuclear cross section describes the interactions of charged leptons(Z) with nuclei 
via virtual photon exchange. The dominant part of the cross section is from nearly real 
photons (negative four- momentum squared Q 2 — ► 0), however, this is the limit where energy 
loss is negligible. Our approach here is to treat the IN — > IX using the deep-inelastic 
scattering formalism, and to use a nucleon structure function F 2 consistent with data over 
the full range of Q 2 . In this approach, both soft physics at low Q 2 and hard perturbative 
physics at high Q 2 are incorporated. The inclusion of the perturbative physics at non-zero 
Q 2 has an effect on the very high energy behavior of f3 nno . 

The standard variables used for l(k)N(p) — > l(k')X scattering include 

q 2 = (k- k' f = -Q 2 (3.1) 

x = (3.2) 
2p ■ q 



p_q 
p ■ k 



(3.3) 



The differential cross section can be written in the form [IE 

da(x,Q 2 ) 4ira 2 F 2 (x,Q 



dQ 2 dx Q 4 x 



1 Mx y a a\ 

1-V—2E (3 - 4) 



/ 2mj\y 2 (l + AM 2 x 2 /Q 2 ) 



Q 2 J 2(1 + R(x,Q 2 )) 
The quantity R is written in terms of Fl and Fi where 



( AM 2 r 2 \ 

F L (x,Q 2 )=il + —^\F 2 (x,Q 2 )-2xF 1 (x,Q 2 ) . (3.6) 

In all of the expressions above, mi is the lepton mass and M is the nucleon mass. We have 
converted the differential cross section to a dependence on y and Q 2 and used the following 
limits of integration: 

QLn <Q 2 < ZMEy - ((M + m^f - M 2 ) (3.7) 
Z/min < V < 1 - mi/E . (3.8) 

where Q 2 min ~ m 2 y 2 /(l - y) and y min ~ ((M + m n ) 2 - M 2 )/(2ME). 

The structure function F L (x, Q 2 ) is proportional to the longitudinal photon- nucleon cross 
section. In the Q 2 — > limit, F L ~ Q A while F 1 —> Q 2 , so R — > 0. We have used R(x, Q 2 ) 
modeled by Badelek, Kwiecihski and Stasto in Ref. |19| for 10 -7 < x < 0.1 and 0.01 
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GeV 2 < Q 2 < 50 GeV 2 . For x > 0.1, the parameterization of Whitlow et al. pO] is used. 



There is no evidence for target dependence of R. The photonuclear /3 nuc for muons evaluated 
with the R parameterized by Badelek et al. differs only by a few percent from that calculated 
with R = 0. Consequently, in what follows, we set R = 0. 

The nuclear structure function depends on the particular target. The attenuation of 
quark density in a nucleus has been observed in deep inelastic lepton scattering (DIS) from 
nuclei at CERN and Fermilab [2~2\ energies in the region of small values of x and Q 2 



The data, taken over a wide kinematic range 10" 5 < x < 0.1 and 0.05 GeV 2 < Q 2 < 100 
GeV 2 , show a systematic reduction of nuclear structure function F^(x,Q 2 )/A with respect 



to the free nucleon structure function 



F 2 N (x,Q' 



a(A, x, 



. We define the shadowing ratio by 

x,Q 2 ) 



AF 2 N \ 



x,Q 2 



(3.9) 



The origin of nuclear shadowing effect is coherent scattering of the virtual photon off 
the nucleons inside the nucleus in the small x region. The coherent multiple scattering 
can be most conveniently handled by the Glauber diffractive approximation model ||23| . 
The interactions of the virtual photons with nucleons is modeled differently for high Q 2 
and low Q 2 . At low Q 2 (< 1 GeV 2 ), calculations of shadowing |24| have often used the 



vector- meson-dominance (VMD) model [ES|, in which the virtual photon interacts with the 
nucleons via its hadronic fluctuations, namely the p, u and mesons. At higher Q 2 , the 
picture is instead that the virtual photon interacts with partonic components of the nucleons 



via its quark- ant iquark pair (qq) color-singlet fluctuation p6 |. Given the high nucleon and 
parton densities, the quarks and gluons that belong to different nucleons in the nucleus will 
recombine and annihilate, leading to the so-called recombination effect first suggested by 
Gribov, Levin and Ryskin [2"7| and later proven by Mueller and Qiu |28| . The net effect of 
either mechanism is that F 2 is lower than one would expect by naive superposition (AF 2 N ) 
and weakly dependent on Q 2 in the range of interest for the photonuclear cross section. We 
take a Q 2 independent function of x and A consistent with the Fermilab E665 data taken 
in the kinematic range of 0.0001 < x < 0.56 and 0.1 < Q 2 < 80 GeV 2 
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A-° A x < 0.0014 , 

a(A, x,Q 2 ) ~ a(A, x) = { A 069 lo ^o -+0.097 .0014 < a; < 0.04 

1 0.04 < x . 



(3.10) 



The structure function F 2 is approximated by 

F 2 A = a(A,x)^(F 2 p + F 2 n 



(3.11) 



a(A,x)A-(l 



P{x))F p 2 



2.45a; 2 — 2.35a; 3 + a; 4 describes the ratio F 2 / F 2 , 



assuming Z = A/2. Here P(x) = 1 — 1.85a; 
parameterized by the BCDMS experiment 

The quantity F 2 (x, Q 2 ) is extracted in a variety of experiments in a range of < Q 2 < 
5000 GeV 2 and 5 x 10~ 6 < x < 1, though kinematic limits restrict the range of Q 2 and 
i in a given experiment. The differential cross section must be integrated from Q 2 = 0, 
where the perturbative QCD description of F 2 is not valid, to values of Q 2 where QCD is 
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valid. Consequently, a parameterization of F% consistent with all the data is most useful 
for our purposes. The parameterization of F% used here is the one by Abramowicz, Levin, 
Levy and Maor (ALLM) [RD|. The ALLM parameterization involves two terms: a pomeron 



contribution and a reggeon contribution. Parameters are used to fit all data available from 
the pre-HERA era as well as HI and ZEUS data published through 1997. The specific form 
with parameters is detailed in the Appendix B. 

Eq. (3.4) shows the alternative to the BB formula for da/dy which appears in Appendix 
A. As an initial comparison of the two approaches, we show in Fig. 1 the cross section for 
real photon-nucleon scattering, as a function of incident photon energy, 

, s , An 2 aF? 

'M-gH,-? 3 - < 3 - 12 > 

indicated by the solid line, and the Bezrukov-Bugaev cross section (dashed line). Photon- 



proton data collected in Ref. |3T] are also shown. Our cross section agrees with the BB 
parameterization at energies below E 1 ~ 10 4 GeV, however the BB cross section increases 
more quickly with energy. At E 1 = 10 9 GeV, our cross section is 0.40 mb, while the BB 
cross section is 0.58 mb. 

For the calculation of the lepton propagation in rock, it is useful to compare the jA cross 
section for standard rock (A = 22). Here we include the nuclear shadowing effects. Our 
results are shown in Fig. 2 with the solid line. The BB cross section is the dashed line. We 
note that our results are in agreement with the BB parameterization for jA, namely 

a(jA) = Aa(-fN)[0.75G(x) + 0.25] , (3.13) 

where x, G(x) and a^N) appear in Eq. (All), over a wide energy. The largest deviation 
occurs at the lowest photon energy shown, where the photo-nuclear contribution is least 
important for charged lepton energy loss. 



IV. MUON ENERGY LOSS, SURVIVAL PROBABILITY, AND RANGE 

The expected average muon energy loss in traversing a material of depth AX is charac- 
terized by (dE/dX), indicated in Eq. (1.1). The results for the standard bremsstrahlung, 
pair production and photonuclear (BB) differential cross sections summarized by LKV, and 
our results of using the ALLM differential cross section are shown in Fig. 3. Our result for 
pnuc fogging t diverge from the standard values at E ~ 10 6 GeV, and is a factor of about 
1.6 higher at E — 10 9 GeV. In terms of the total (3, the total with ALLM contributions is a 
factor of 1.15 larger than with the BB photonuclear contribution at E — 10 9 GeV. 

To explore the effect of the slightly larger value of f3 nnc for muons, we have evaluated the 
muon survival probability in standard rock. The survival probability P(E, X) for a muon to 
survive to a depth X given incident energy E incorporates the effects of fluctuations due to 
radiation. In Fig. 4, we show our muon survival probabilities (solid lines) for E = 10 3 — 10 9 
GeV, in decades of energy, versus survival depth X (in km.w.e.) for standard rock (A = 22 
and p = 2.65 g/cm 3 ). We have taken E min = 1 GeV in the Monte Carlo. Using the 
LKV defaults in our Monte Carlo program yields the dashed lines, which agree with the 



Lipari-Stanev result in Ref. (17 ] 
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The two sets of survival probabilities translate to average muon ranges with incident 
energy E and final energy E min = 1 GeV. The average range (R(E)) is defined by 

(R(E))= / dXP(E,X) . (4.1) 
Jo 

The average ranges for our calculation are shown in Fig. 5 by the solid lines. The standard 
LKV ranges that we have calculated using the same muon transport Monte Carlo simulation 
are shown with dashed lines. At E — 10 9 GeV, the two calculations differ by only 5%. The 
deviation from the standard calculation increases with energy. 

V. TAU ENERGY LOSS, SURVIVAL PROBABILITY, AND RANGE 

Essentially the same procedure for calculating muon energy loss can be applied to the 
tau lepton, with the important modification that the tau has a decay length considerably 
shorter than the muon decay length. We have evaluated the pair production, bremsstrahlung 
and ionization energy loss according to the formulas in Appendix A for tau leptons, and we 
have evaluated the tau photonuclear differential cross section using the ALLM expression 
in Eq. (3.4). Fig. 6 shows the tau energy loss contributions to (3 for standard rock. In this 
figure, f3 is plotted on a logarithmic scale because the bremsstrahlung is very suppressed 
relative to the other contributions, due to the much heavier lepton mass. The photonuclear 
contribution to f3 dominates above E ~ 10 5 GeV. 

In the absence of energy loss, the tau survival probability corresponding to the curves in 
Fig. 4 are exponentials of the form: 

(5.1) 

where 7 = E/m T is the Lorentz gamma factor, cr = 86.93 /iin is the tau decay length and 
p is the material density. The average range with no energy loss, as defined by Eq. (4.1), is 
just 

(R(E)} decay = ICTp . (5.2) 

For incident energies between 10 3 GeV and 10 9 GeV, we have made a Monte Carlo 
evaluation of tau energy loss in water, rock and iron including the electromagnetic energy 
loss mechanisms as described in the Appendix. 

At low energies, the decay length of the tau is short and energy loss is relatively unimpor- 
tant, so the survival probability is just the exponential in Eq. (5.1). The survival probability 
is increased at fixed depth as the energy of the tau increases. The survival probability curves 
can be put on the same plot by using P(E, X) versus X/E as shown in Fig. 7 for water. In 
Fig. 7, the decay distribution Eq. (5.1) is indicated by a dashed line. A solid line overlays it, 
which is the probability distribution, as computed by our Monte Carlo, with incoming tau 
energy E = 10 3 GeV. The lower solid line is the survival probability for incident tau energy 
E = 10 9 GeV. The probabilities are evaluated with a minimum tau energy of 50 GeV. By 
changing the minimum energy to 100 GeV, we see no change in the probability distributions 
for incident tau energies E = 10 3 — 10 9 GeV. 



P{E,X) =exp 



X 

jcrp 
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In Fig. 8, we compare the tau decay length with the range of tau with incident energy 
E and final energy E min = 50 GeV. The dashed line shows the tau decay length (Eq. (5.2)), 
while the solid line shows the evaluation of Eq. (4.1) including electromagnetic energy loss 
for propagation in water. The deviation from the simple gamma factor scaled decay length 
starts at about 10 8 GeV in water, and by 10 9 GeV, the average range is 35% lower. The 
dot-dashed and dotted lines in Fig. 8 show the tau ranges in standard rock and iron. For 
tau propagation in iron, the tau range is an order of magnitude shorter than its decay length 
for an incident energy of 10 9 GeV. 

Our results differ from the estimate of the range by Fargion in Ref . , mainly because 
of the inclusion of the photonuclear and pair production processes in electromagnetic energy 
loss. The estimate in Ref. relies on an approximate solution to Eq. (1.1), where a and 
(3 are assumed to be energy independent, and (3 is approximated numerically by rescaling 
the muon pair production (3. The resulting (3 is more than an order of magnitude smaller 
than the (3 values we obtain, for example, in Fig. 6 for rock. As a consequence, whereas in 
Ref. ||32|| , electromagnetic energy loss is never relevant, we have shown that accounting for 



photonuclear, bremsstrahlung and pair production energy loss mechanisms reduces the tau 
range in water beginning at E ~ 10 8 GeV, at even lower energies for more dense materials. 



VI. DISCUSSION 

We have re-evaluated the muon energy loss due to photonuclear interactions using the 
recent HERA results for the real and virtual photon-nucleon scattering. We have used the 
ALLM parameterization of the electromagnetic structure function F 2 (x, Q 2 ) to evaluate the 
photonuclear cross section including the Q 2 3> region. As compared with the previous 
Bezrukov and Bugaev result, our approach yields a change in (3 for photonuclear interactions. 

In our evaluation of energy loss and effective ranges, except for tau decays, we have 
ignored the effects of weak interactions on the propagation of muons and taus. To set the 
scale for charged lepton energy loss via weak interactions, 

flweak - N A (y}ai N (E) . 

For (y) ~ 0.2 and using <r lN ~ a uN from Ref. ||, we find that for j3 wea k to be larger than 
a typical high energy electromagnetic (3 ~ 10~ 6 cm 2 /g, the charged lepton energy must be 
larger than ~ 10 16 GeV. Below that energy, weak interactions do not play a significant role 
in charged lepton energy loss. 

Apart from energy loss, weak interactions also contribute to charged lepton disappear- 
ance through charged-current interactions. The tau charged-current interaction length is 
comparable to its decay length at energies above 10 10 GeV in water. We have shown in Fig. 
8 that already at 10 8 GeV, taus lose energy in water due to electromagnetic interactions. 
For tau energies above 10 10 GeV in water, it will be an interplay of electromagnetic energy 
loss, decay and charged-current weak interaction disappearance of taus that will dictate tau 
effective ranges. Even in lead, which has a density of 11.35 g/cm 3 , the decay length is shorter 
than the charged-current interaction length for energies below a few times 10 9 GeV, so our 
neglect of weak interactions below E = 10 9 GeV is a good approximation. 
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The results presented here are not significantly modified by the Landau- Pomeranchuk- 
Migdal (LPM) effect [33|. The LPM effect arises from low-momentum transfer lepton scat- 



tering, where bremsstrahlung should, in principle, be evaluated including coherence effects 
due to scattering in the medium. The effect has been measured in the case of electron 
scattering [|33J and summarized in Ref. [j3"5| , but is suppressed for heavier leptons. The 



suppression factor can be written crudely in terms of the photon energy q — y ■ E as [35 



when the argument of the square root is positive (y < E/(E + E LPM )), where 

E LPM = 1.38 x 10 13 ^— • X . (6.2) 
cm 

The radiation length Xo is 36.1 cm for water and 10.0 cm for rock. For the lepton energies 
of interest here, E < 10 9 GeV, this means that LPM suppression for muons occurs for 
y < 7 x 10~ 5 (cm/Xo). In practice, this affects the continuous "soft" term in Eq. (2.1). The 
approximation for S(y, E) leads to a decrease in the soft term by at most 0.1% at the highest 
energies considered here. This is consistent with, for example, the more detailed evaluation 
of /3 hrem by Polityko et al. in Ref. [36|, where deviations from standard f3 hrera evaluations 



occur above E = 10 11 GeV. The LPM effect is even further suppressed for taus, since Elpm 
scales like the lepton mass to the fourth power. 

The change to (3 from using the ALLM photonuclear interaction cross section increases 
with energy and reaches 60% for E = 10 9 GeV for muons. Still, for muons, the photonu- 
clear processes do not dominate (3 for E < 10 9 GeV, and the overall effect on dE/dX is 
smaller. The muon survival probabilities and average muon ranges are obtained from the 
one dimensional Monte Carlo approximation. We find our results to be within 5% from the 
standard LKV value with incident energy E < 10 9 GeV and final energy E min = 1 GeV. 

We have evaluated the tau energy loss including the photonuclear processes as well as the 
ionization, pair production and bremsstrahlung energy loss. At energies above ~ 10 8 GeV, 
tau interactions become important in water and the range becomes significantly shorter. 
The effect appears at lower energies for more dense materials, for example, at tau energies 
less than 10 7 GeV for iron. We thus expect a decrease in the observed tau flux and tau 
neutrino energy relative to that expected if these effects had be ignored. The effect is most 
important for high energy taus arriving from directions just below the horizon. For directions 
subtending a significant portion of the Earth, energies will be degraded to less than 10 6 GeV 
by repeated neutrino conversion to tau and tau decay. The last few steps in that process will 
occur at energies where dE/dX is not important. Near the horizon, the interaction length 
for charged current conversion becomes comparable to the length of the chord through the 
Earth, and the produced taus may arrive directly in the detector depending on their range. 
For E > 10 s GeV, the tau range is decreased in water and such taus will either not arrive or 
arrive at lower energies. Correspondingly, for taus produced in a detector, measured dE/dX 
would be higher than expected without the revised photonuclear effects calculated here. 
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APPENDIX A: ENERGY LOSS FORMULAS 



For completeness, we reproduce the standard formulas for ionization, pair production and 
bremsstrahlung energy loss. These are also collected in the unpublished work of Lohmann, 
Kopp and Voss ||. The Bezrukov and Bugaev (BB) parameterization of photonuclear energy 
loss is also included below. 

The constants are defined as follows: 

a= 1/137 (Al) 
N A = 6.023 x 10 23 , Avogadro's number 

Z— atomic number, A = atomic weight of medium 
m e , mi rest masses of electron and muon or tau 
A e = 3.8616 x 1CT 11 cm, Compton wavelength of the electron 
e= 2.718 



To account for energy loss on atomic electrons, the factor Z 2 in the bremsstrahlung and 
pair production formulas below should be replaced by Z(Z + 1). 



1. Ionization Energy Loss 

The ionization loss for a incident muon or tau of energy E and momentum p is given by 



the Bethe-Bloch formula 11 



where 



— = a 2KN A X e -^ (In /2(z) - 2/3 + — - 5(X) j (A2) 



/3 = p/E, 
7 = E/mi 



The quantity 



E' = 2m e — 1 (A3) 



is the maximum energy that can be transfered to the electron. A density correction is 
parameterized by 5(X) for X = log l0 (/3 , ~f) by Sternheimer et al. JT2| as 



S(X)= 4.6052X + a(X l - X) m + C X < X < X 1 (A4) 

5(X)= 4.6052A + C X > X 1 . 



The values of X , X\, a, m, C and the mean ionization potential I{Z), for rock and water 
are given in Table 2. 
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2. Bremsstrahlung Energy Loss 



Following Petrukhin and Shestakov's evaluation of the form factors in the Bethe-Heitler 
bremsstrahlung differential cross section |T3|], the differential cross section used here is: 



where 



dy \ mi J y\3 3 J 



0(5)= In 
<P{5)= In 
<) 



189m; V~ 1/3 
rrt e 

1 + 1893^^-1/3 

- 2 189m; 7~ 2 / 3 
3 m e 

1 + 189v^ 5Z -l7i 

m 2 y 



2E(1 - y) 



The range of y integration is 



Z < 10 



Z > 10 



0<,<l-^v^ 3 . 



(A5) 



(A6) 



(A7) 



3. Pair Production Energy Loss 

Pair production energy loss depends on the asymmetry parameter p given as 



E+-E- , , / 



6mf 



E 2 (l~y), 



'I- 



4m e 
~Ey~ 



(A8) 



given E + and E , the energies of the positron and electron pair. The differential cross 
section is parameterized by pi 



(Pa 



on 



dydp '6n ' y 
where the functions e and (pi are as follows: 



+ —<Pl 

mf , 



(2 + p*)(l +0) + ^(3 + p 2 ) J ln(l + -) + — ^ 
(1 + P 2 )(l + - ~(1 + 2/3)(l - p 2 ) ]ln(l + 



- (3 + P 2 ) 



C(l-p 2 -/?) 
1+^ 



+ (l + 2/3)(l-p 2 



RZ~^(1 + 0(1 + ^e) 1 
±j e — In 1 — — In 

i , 2m eN /SflZ~?(l+g)(l+r e ) 2 

1 + Sy(i-P 2 ) 



i + (|-^) 2 (i + 0(i + K 

2 m ; 



(A9) 



(A10) 
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U= In 

Y e -- 
Y r - 



3 rrif. 



1 + 



2-m ev ^fiZ-?(l+g)(l+y;) 

Ey(l- P 2 ) 

5 - p 2 + 4/3(1 + p 2 ) 



2(1 + 3(3) ln(3 + i) - p 2 - 2/3(2 - p 2 ) 

4 + p 2 + 3/3(l + p 2 ) 
(l + p 2 )(| + 2/3)ln(3 + + l-|p 2 



mil/ 



1 



2(1-2/)' V 2 «W (i-y) * 

The value of i? is R = 189. The range of y integration is 



4m e 



<y<i 



3m; 
4£ 



1/3 



(All) 



4. Bezrukov and Bugaev Photonuclear Energy Loss 



The parameterization of the Bezrukov and Bugaev photonuclear energy loss is |10 

Km\ 2mf 



da 
dy 



a 
27 



A(T jN y 



2mf 



mt 



+ 



m\ l 3 



where 



..... , . ... s mi lm . ,. t 

H - -G{x)^r-± 1 2 

2t U V mj + t At 



+ ±^ln(l + JL) 

m^ 



G(x) 



- 1 + e~ 1 Yl + x) 



x" 



X- 



0.00282^3 ((j 7 jv(E/GeV)/pb) 



a 7iY (^)= H4.3 + 1.6471n 2 (0.0213£;/GeV) pb 
t= 



K= 1 



2 2 

mfy 
2 2 



m-i 



0.54 GeV 2 



m 2 = 1.8 GeV z 



(A12) 



(A13) 



APPENDIX B: ALLM PARAMETERIZATION OF F 2 

The ALLM phenomenological parameterization of F%{x, Q 2 ) has a pomeron and reggeon 
contribution [RDl: 
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Ff (*, g 2 



Q 2 



Q 2 + m 2 

where the two terms have the form 



F?(x,Q 2 )+F«(x,Q 



F?(x } Q 2 ) = c v (t)x a p v{t \l-x) b ^ t) 
F^(x,Q 2 ) = c n (t)x a ^ t \l-x)^ 

The functions c-p(t), etc., are parameterized in terms of 

']n(Q 2 + Q 2 )/A 2 \ 



t = In 

The form for f = c-jz and a-jz is 



hiQg/A 2 



f = fi + f2t h 



for g = op and ap, 



g = gi + {gi- 92) 



1 + #3 



- 1 



and finally, h = b-ji and bp, 



In the expressions for F% and F^ 1 , 



h = h\ + h 2 2 t h3 



Xp 

x-r. 



Q 2 + m 2 P 
Q 2 + m 2 v + W 2 - M 2 
Q 2 + m| 



Q 2 + m\ + W 2 — M 2 ' 
The parameters satisfying the fits described in Ref. [|30| are shown in Table II. 



(Bl) 



(B2) 
(B3) 



(B4) 

(B5) 

(B6) 

(B7) 

(B8) 
(B9) 
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TABLES 



TABLE I. Parameter values in ionization energy loss [12]. 



Material 


J/eV 


-C 


X 


Xi 


a 


m 


Z 


A 


p/(g/cm 3 ) 


Water 


75.0 


3.502 


0.240 


2.800 


0.091 


3.477 


6.6 


11.89 


1.00 


Standard Rock 


136.4 


3.774 


0.049 


3.055 


0.083 


3.412 


11 


22 


2.65 


Iron 


286.0 


4.291 


-0.0012 


3.153 


0.147 


2.963 


26 


55.84 


7.87 



TABLE II. Parameter values for ALLM parameterization of F^{x, Q 2 ). 



C-pi 


0.28067 


cm 


0.80107 


Cp2 


0.22291 




0.97307 




2.1979 


cm 


3.4942 


avi 


-0.0808 


am 


0.58400 




-0.44812 


am 


0.37888 




1.1709 


am 


2.6063 


bpi 


0.60243 


bm 


0.10711 




1.3754 


bm 


1.9386 




1.8439 


bm 


0.49338 


m 2 (GeV 2 ) 


0.31985 


m 2 P {Ge\ 2 ) 


49.457 


m^(GeV 2 ) 


0.15052 


Q 2 (GeV 2 ) 


0.46017 


A 2 (GeV 2 ) 


0.06527 
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FIG. 1. The photon-nucleon cross section as a function of incident photon energy for the BB 
(dashed) and ALLM (solid) parameterizations. Also shown are photon-proton data collected in 
Ref. [30]. 
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FIG. 2. The photon-nucleus cross section standard rock (^4 = 22), as a function of incident 
photon energy E using the ALLM parameterization and Eqs. (3.9) and (3.10) for the shadow 
factor and conversion to nucleon structure function. 
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FIG. 3. The f3 value for muon in standard rock (A = 22), including bremsstrahlung (solid line), 
pair production (dashed) and photonuclear (dotted) interactions. 
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FIG. 4. Muon survival probabilities in rock using BB differential cross section for the photonu- 
clear term (dashed) and the ALLM differential cross section (solid). 
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FIG. 5. Average muon range in standard rock (km.w.e. depth), for incident muon energy E, 
final muon energy Ef > 1 GeV, using the standard LKV treatment of energy loss, including the 
BB differential cross section (dashed) and substituting the ALLM photonuclear calculation (solid). 
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FIG. 6. The (3 value for tau in standard rock (A = 22), including bremsstrahlung (solid line), 
pair production (dashed) and photonuclear (ALLM) (dotted) interactions. 
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X/pE [cm/GeV] 

FIG. 7. Decay distribution(dashed line) and survival probability curve for incoming tau energy 
E = 10 3 GeV (upper solid line) and E = 10 9 GeV (lower solid line) in water (p = 1 g/cm 3 ). 
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FIG. 8. Tau decay length (dashed line) and the average tau range in water (solid line), in rock 
(dot-dashed line) and in iron (dotted line), for incident tau energy E, final tau energy larger than 
Emin = 50 GeV including electromagnetic energy loss. 
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